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One possibility for Beyond Standard Model physics is a new strongly-interacting gauge theory. 
One way to determine if a non-abelian gauge theory is QCD-like or conformal is to measure the 
running of the renormalized gauge coupling. We define the renormalized coupling from Wil- 
son loop ratios, and measure these ratios via lattice simulations. We test this method in SU(3) 
pure gauge theory and show some first results for simulations with dynamical fermions in the 
fundamental representation. 
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1. Running coupling 

There is great interest in the possibility that Beyond Standard Model physics might take the 
form of new strongly coupled gauge theories fll|]-[|2ll], one example being technicolor. For model 
building, it is necessary to distinguish conformal theories from those with a mass gap like QCD. 



There have been a number of lattice studies of this issue [|22|]-[|33[]. One signal is the running 
coupling, which has an infrared fixed point in a conformal theory. We developed our method 
for the running coupling and its beta function from the definition of the continuum renormalized 
coupling constant using the second derivative of R x T Wilson loops p4|]. This coupling runs with 



the size of the Wilson loop in infinite volume Q34|]. We generalized this definition to a finite volume 
L keeping R/L fixed and run the coupling with L as in the Schrodinger functional method. A 



similar method was developed independently in Q33| ] 



Consider Wilson loops W(R,T,L), where R and T are the space-like and time-like extents 
of the loop, and the lattice volume is L 4 (all dimensionful quantities are expressed in units of the 
lattice spacing a). A renormalized coupling can be defined by 

^ R i L ^--mTL)^sf in{v,{RJ - L)) ^ <«> 

where for convenience the definition will be restricted to Wilson loops with T =R, and (...) is the 
expectation value of some quantity over the full path integral. This definition can be motivated 
by both renormalized and bare perturbation theory, where the leading term is the tree-level cou- 
pling. The renormalization scheme is defined by holding R/L to some fixed value. The quantity 
k(R/L) can be calculated from Wilson loop expectation values using perturbation theory. This is 
done numerically on finite lattices, hence k contains lattice artifacts which vanish as L — > oo. The 
role of lattice simulations is to measure the expectation values non-perturbatively. On the lattice, 
derivatives are replaced by finite differences, so the renormalized coupling is defined to be 

g 2 ((R + l/2)/L,L) = -±-^R + l/2) 2 X (R+l/2,L) , 



X(R + l/2,L) = -\n 



W{R + 1,T + 1,L)W(R,T,L) 
W(R + l,T,L)W(R,T + l,L) 



\T=R J 



(1.2) 



where % is the Creutz ratio []35[], and the renormalization scheme is defined by holding the value of 
r=(R + l/2)/L fixed. 

With this definition, the renormalized coupling g 2 is a function of the lattice size L and the 
fixed value of r. The continuum limit corresponds to L — > oo, where the physical length scale L p h ys 
is held fixed while the lattice spacing a — > 0. The coupling is non-perturbatively defined, as the 
expectation values are calculated via lattice simulations, which integrate over the full phase space 
of the theory. One starts the RG flow from some reference physical point L p h ys ,o which is set by the 
choice e.g. g (r,Lphys,o) = 0.8. m a QCD-like theory, g 2 increases with increasing L p h ys flowing in 
the infrared direction. In a conformal theory, g 2 flows towards some non-trivial infrared fixed point 
g as Lphy S increases, whereas in a trivial theory, g 2 decreases with L p h ys . 

One way to measure the running of the renormalized coupling in the continuum limit is via 
step-scaling. The bare lattice coupling is defined in the usual way /3 = 6/gg as it appears in the 
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Figure 1: (Left) The rescaled Creutz ratio on a 28 4 lattice at j3 = 6.99. We interpolate the data linearly to 
r = (R + l/2)/L = 0.25, giving a chi squared per degree of freedom 2.0/4. (Right) The measured coupling 
g 2 {2Li,fii) for 2Li = 20,24,28 and 32, where J3; is tuned such that g 2 (Z,,-,j3,-) = 1.44. A linear continuum 
extrapolation gives g 2 (2L phys ) = 1.636(23), with £ 2 /dof = 0.57/2. 

lattice action. On a sequence of lattice sizes L\,L2, ...,L n , the bare coupling is tuned on each lattice 
so that exactly the same value g 2 (r,L,-, /$,) = g 2 (r,Lphys) is measured via simulations. Next a new 
set of simulations is performed, on a sequence of lattice sizes 2Li,2L}, ...,2L n , using the corre- 
sponding tuned couplings j3i, /fe, From the simulations, one measures g 2 (r,2L,-, /3,), which 
vary with the bare coupling i.e. the lattice spacing. These data are extrapolated to the contin- 
uum as a function of 1/L?. This gives one blocking step g 2 (r,L p h ys ) — ► g 2 (r,2L p h ys ) in the contin- 
uum RG flow. The whole procedure is then iterated. The chain of measurements gives the flow 
g 2 (r,L phys ) -> g 2 (r,2L phys ) -> g 2 (r,4L phys ) -> g 2 (r,8L phys ) -> as far as is feasible. One is free to 
choose a different blocking factor, say L p h ys — > (3/2)L p h ys , in which case more blocking steps are 
required to cover the same energy range. 

2. St/ (3) pure gauge theory 

As our first test of this method, we study SU (3) pure gauge theory in four dimensions. We 
simulate using the standard Wilson lattice gauge action, with a mixture of five over-relaxation 
updates for every heatbath update. We define the renormalization scheme with the fixed value 
r = 0.25, for brevity we omit the label r in the renormalized coupling. In this R/L range, the Creutz 
ratio can be accurately measured, and the geometric factor k converges quickly to its continuum- 
limit value. In the pure gauge theory test, we actually use the finite L values of k, as this may 
remove some of the cutoff dependence of the renormalized coupling. For the RG flow we choose 
the blocking step L — > 2L. We simulate on small lattices of size L = 10, 12, 14, 16, 18,20 and 22, and 
the corresponding doubled lattices 2L = 20,24,28,32,36,40 and 44. Wilson loop measurements 
are separated by 1000 sweeps, which we find is sufficient to generate statistically independent 
configurations. 

To tune /3, for each small lattice size L,, we typically run separate simulations at 5 - 10 different 
/3 values in the relevant range of renormalized coupling g 2 . Each of these runs contains 300 - 500 
measurements i.e. up to 5 x 10 5 sweeps. Simulating on the doubled lattices at the tuned j8, values, 
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Figure 2: (Left) The continuum extrapolations of four discrete RG steps. (Right) The RG flow g 2 (L p h ys ), 
combining analytic lattice perturbation theory and the simulation results. The running starts at g 2 (L p hy S .o) = 
0.825. There is excellent agreement with continuum 2-loop running, at the strongest coupling, the simulation 
results begin to break away from perturbation theory. 



we generate between 200 and 1000 measurements each, typically more than 500. The signal of the 
Creutz ratio disappears into the noise as the size of the Wilson loop increases. One way to suppress 
the noise in measurements is to gauge fix the configurations to Coulomb gauge, and replace the 
thin-link Wilson loop with the correlator of the products of the time-like gauge links [3^, 37]. Note 
that gauge fixing is not implemented in the actual Monte Carlo updating algorithm. An alternative 
method to suppress noise is to smear the gauge links and measure the fat-link Wilson loop operator. 
In the pure gauge theory test, we use the gauge-fixing method, in the dynamical fermion simulations 
we describe later, we use the smearing method. These improvement methods do not correspond to 
calculating the original thin-link Wilson loop operator. 

We show in Fig. [I] (left) a typical result for the rescaled Creutz ratio {R + 1 /2) 2 X- The doubled 
lattice is 28 4 and the bare coupling j3 = 6.99 is tuned from simulations on 14 volumes. Errorbars 
are calculated using the jackknife method. The renormalized coupling is defined at the point r = 
(/?+ l/2)/L = 0.25, corresponding to (R + 1/2) = 7 at this lattice size. We interpolate the data 
linearly to this point, obtaining # 2 /dof= 2.0/4. The data at different R are highly correlated, being 
measured on the same gauge configurations. To calculate an error, we bin the gauge configurations, 
and analyze and interpolate separately each bin, giving a distribution of interpolated values. An 
example of the step-scaling method is shown in Fig. [I] (right). The bare couplings are tuned such 
that g 2 (L ; ,j8,) = 1.44 for U = 10,12,14 and 16, the figure shows the data g 2 (2L/,j8,) for 2L,- = 
20,24,28 and 32. The leading lattice artifacts in the Creutz ratio are expected to be of order 
&{a 2 ). This corresponds to ff(\/L 2 ), since the physical lattice size is La. Extrapolating linearly 
in 1/ (2L) 2 gives g 2 (2L p h ys ) = 1.636(23) and # 2 /dof = 0.57/2. For a systematic error, we omit one 
data point at a time and repeat the extrapolation. Combined in quadrature with the statistical error, 
our continuum result is g 2 (2L phys ) = 1.636(25). 

We iterate the procedure, giving four discrete RG steps as shown in Fig. || (left). At stronger 
coupling, we need lattices up to 2L, = 44 for the continuum extrapolation. The use of the finite 
L values of k does not appear to reduce the cutoff effects. The continuum RG flow is shown in 
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Figure 3: The renormalized coupling g 2 (L,fi) at fixed bare coupling for Nf = 16 fundamental flavors. For 
g 2 (L : fi) > 0.5 the renormalized coupling decreases with increasing L, for g 2 (L,fi) < 0.5 the renormalized 
coupling is independent of L within the errors. This is consistent with the existence of an infrared fixed 
point. 

Fig. U (right). At weak coupling, we use analytic/numeric lattice perturbation theory to calculate 



the Wilson loop ratios in finite volumes [38, 39]. The Wilson loops are calculated to 1-loop in 



the bare coupling in finite volume. The series is reexpanded in the boosted coupling constant 
at the relevant scale of the Creutz ratio. Step-scaling of the finite volume ratios can be used in 
exactly the same way as for the simulations results, to determine the RG flow in the continuum. 
The analytic RG flow starts at the reference point g 2 (Lphys.o) = 0.825. At weak coupling there is 
complete agreement with 2-loop perturbation theory. We connect lattice perturbation theory to the 
simulation results by matching the flows at g 2 (L p h ys ) = 1.44, where the simulation RG flow begins. 
There is continued agreement with 2-loop perturbation theory at even stronger coupling, only at 
the strongest coupling do we see deviation from the perturbative flow. 

3. Fundamental fermions 

We next study SU (3) gauge theory with fermions in the fundamental representation. For 
Nf = 16 flavors, 2-loop perturbation theory predicts the theory is conformal with an infrared fixed 
point g* 2 0.5. Because of the computational expense of step-scaling with dynamical fermions, 
in this pilot study we have not yet extrapolated to the continuum limit. The running coupling 
therefore is still contaminated with finite cutoff effects. If the linear lattice size L is large enough, 
the trend from the volume dependence of g 2 (L, /3 ) should indicate the location of the fixed point. 
For g 2 (L,f}) > g* 2 we expect decrease in the running coupling as L grows, although the cutoff of 
the flow cannot be removed above the fixed point. Below the fixed point with g 2 (L,fi) < g* 2 we 
expect the running coupling to grow as L increases and the continuum limit of the flow could be 
determined. The first results are shown in Fig. |3[ We use stout-smeared [ 40 ] staggered fermions 



lj, |42J and the RHMC algorithm, simulating at quark mass m q = 0.01, with some runs at m q 
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0.001 to test that the mass dependence is negligible. For the Wilson loop ratios, we smear the gauge 
fields and measure the fat-link Wilson operator. Our experience in the pure gauge theory test is that 
cutoff dependence is not reduced using the finite L value of k, hence we use the infinite volume 
k value to convert the Wilson loop ratios to a renormalized coupling. The results are consistent 
with the above picture. For g 2 (L,/3) > 0.5, the cutoff dependent renormalized coupling decreases 
with L. For g 2 (L,j8) < 0.5, the renormalized coupling is L-independent within errors. The theory 
appears conformal but precise determination of the conformal fixed point requires further studies. 
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